O 



c 



in 

00 



o 



Are motorways rational from slime mould's point of view? 

Andrew Adamatzky 1 , Selim Akl 2 , Ramon Alonso-Sanz 3 , Wesley van Dessel 4 , Zuwairie Ibrahim 5 , 
Andrew Ilachinski 6 , Jeff Jones 1 , Anne V. D. M. Kayem 7 , Genaro J. Martinez 8 , Pedro de Oliveira 9 , 
Mikhail Prokopenko 10 , Theresa Schubert 11 , Peter Sloot 12 , Emanuele Strano 13 , and Xin-She Yang 14 

1 University of the West of England, Bristol, UK 

2 Queen's University, Kingston, Ontario, Canada 



4 



3 Universidad Politecnica de Madrid, Madrid, Spain 
Scientific Institute of Public Health, Brussels, Belgium 
Universiti Teknologi Malaysia, Johor Darul Takzim, Malaysia 
6 Center for Naval Analysis, Alexandria, USA 

7 University of Cape Town, Cape Town, South Africa 

8 Universidad Nacional Autonoma de Mexico, Mexico 
^^ 9 Mackenzie Presbyterian University, Sao Paulo, Brazil 

10 CSIRO Information and Communication Technologies Centre, Sydney, Australia 
^ n Bauhaus-Universitat Weimar, Weimar, Germany 

12 University of Amsterdam, Amsterdam, The Netherlands 
^i 13 eppl^ Lausanne, Switzerland 

£f*} 14 National Physical Laboratory, Teddington, United Kingdom 



Abstract. We analyse the results of our experimental laboratory approximation of motorways net- 
works with slime mould Physarum polycephalum. Motorway networks of fourteen geographical areas 
are considered: Australia, Africa, Belgium, Brazil, Canada, China, Germany, Iberia, Italy, Malaysia, 
Mexico, The Netherlands, UK, USA. For each geographical entity we represented major urban areas by 
oat flakes and inoculated the slime mould in a capital. After slime mould spanned all urban areas with 
a network of its protoplasmic tubes we extracted a generalised Physarum graph from the network and 
compared the graphs with an abstract motorway graph using most common measures. The measures 
employed are the number of independent cycles, cohesion, shortest paths lengths, diameter, the Harary 
index and the Randic index. We obtained a series of intriguing results, and found that the slime mould 
approximates best of all the motorway graphs of Belgium, Canada and China, and that for all enti- 
ties studied the best match between Physarum and motorway graphs is detected by the Randic index 
(molecular branching index). 
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1 Introduction 



The increase of long-distance travel and subsequent reconfiguration of vehicular and social networks [35] 
requires novel and unconventional approaches towards analysis of dynamical processes in complex transport 
networks [15] , routing and localisation of vehicular networks [43] , optimisation of interactions between differ- 
ent parts of a transport network during scheduling of the road expansion and maintenance [53] , and shaping 
of transport network structure [17]. "The concept of a network is useful, but it may be misleading, because 
it gives the impression that it arises from a coherent planning process, whereas in practice almost all road 
networks have evolved gradually, being added to and altered over time to meet the ever-changing needs and 
demands of travellers. (...) even our more recently built network of roads, the motorways designated for 
exclusive use of motor vehicles, has been result of nearly a century of thinking and rethinking about what 
should be built." [21] 

Modern motorway networks are based on the millennium- long emergence of roads. First there were pre- 
historic trackways — mesolithic footways established by early men 'seeking most direct and convenient 
alternatives by process of trial and error' |16)21j or based on the routes selected by the animals [54 , tim- 
ber trackways and droveways (used primarily by cattle). Further development of roads was country-specific. 
Talking about England we can speculate that Romans were building roads along pre-existing trackways and 
possibly along Ridgeways. In the 1700s turnpikes were established. They were based on pre-existing roads 



with few local diversions [18] . The turnpikes were then substituted by single carriageways, dual carriageways 
and, finally, motorways [21] . Thus backtracking the history of motorways we arrive to pathways developed 
by living creatures. How were the pathways developed? Were they efficient? More or less accurate answers 
could be found by imitating the road network development with living substrates. 

While choosing a biological object to imitate the growth of road networks we want it to be experimental 
laboratory friendly, easy to cultivate and handle, and convenient to analyse. Ants would indeed be the first 
candidate. They do develop trails very similarly to pre-historic people. The great deal of impressive results has 
been published on ant-colony inspired computing |23|49j . However ant colonies require substantial laboratory 
resources, experience and time in handling them. Actually very few, if any, papers were published on exper- 
imental laboratory implementation of ant-based optimisation (also there is the issue of what happens if ant 
collide on the paths, or fight), the prevalent majority of publications being theoretical. There is however an ob- 
ject which is extremely easy to cultivate and handle, and which exhibits remarkably good foraging behaviour 
and development of intra- cellular transport networks. This is the Plasmodium of Physarum polycephalum. 
Plasmodium is a vegetative stage of acellular slime mould P. polycephalum, a single cell with many nuclei, 
which feeds on microscopic particles [50]. When foraging for its food the Plasmodium propagates towards 
sources of food particles and microbes, surrounds them, secretes enzymes and digests the food. Typically, the 
Plasmodium forms a congregation of protoplasm covering the food source. When several sources of nutrients 
are scattered in the Plasmodium's range, the Plasmodium forms a network of protoplasmic tubes connecting 
the masses of protoplasm at the food sources. A structure of the protoplasmic networks is apparently optimal, 
in a sense that it covers all sources of nutrients and provides a robust and speedy transportation of nutrients 
and metabolites in the Plasmodium's body [38 39|40j . 

Motorway networks are designed for efficient vehicular transportation of goods and passengers, protoplas- 
mic networks are developed for efficient intra-cellular transportation of nutrients and metabolites. Is there a 
similarity between these two networks? 

To uncover analogies between biological and human-made transport networks and to project behavioural 
traits of biological networks onto development of vehicular transport networks we conducted a series of 
experimental laboratory studies on evaluation and approximation of motorway networks by P. polycephalum 
in fourteen geographical regions: Africa, Australia, Belgium, Brazil, Canada, China, Germany, Iberia, Italy, 
Malaysia, Mexico, The Netherlands, UK, and USA [2] -[14], [51]. We represented each region with an agar 
plate, imitated major urban areas with oat flakes, inoculated Plasmodium of P. polycephalum in a capital, and 
analysed structures of protoplasmic networks developed. For all regions studied in laboratory experiments [2]- 
[14], [51], we found that, the network of protoplasmic tubes grown by Plasmodium matches, at least partly, 
the network of human-made transport arteries. The shape of a country and the exact spatial distribution of 
urban areas, represented by sources of nutrients, may play a key role in determining the exact structure of 
the Plasmodium network. In the present paper we aim to answer two principal questions. What measures, 
apart from straightforward comparison of edges between motorway and Plasmodium networks, are reliable 
indicators of matching? Which regions have the most 'Physarum friendly' motorway networks, i.e. show 
highest degree of matching between motorways and protoplasmic networks along several measures? In the 
course of investigation we got the answers to these questions and also obtained a few quite intriguing results 
on hierarchies of motorway and protoplasmic networks based on measures and topological indices. 

2 Experimental 

In laboratory experiments [2]-[l4],[5T] we considered fourteen regions: Australia, Africa, Belgium, Brazil, 
Canada, China, Germany, Iberia, Italy, Malaysia, Mexico, The Netherlands, UK, and USA. Agar plates, 2% 
agar gel (Select agar, Sigma Aldrich), were cut in a shape of any particular region and placed in polyestyrene 
square Petri dishes 120 x 120 mm or 220 x 220 mm. For each region we choose the most populated urban 
areas U scaled down locations of which were projected onto agar gel. Numbers of the areas selected for each 
country are as follows: Africa n = 35, Australia n = 25, Belgium n = 21, Brazil n = 21, Canada n = 16, 
China n — 31, Germany n — 21, Iberia n = 23, Italy n — 11, Malaysia n = 20, Mexico n — 19, Netherlands 
n — 21, UK n = 10, USA n = 20, see detailed configurations in [2]-[T4],[5T]. 

At the beginning of each experiment a piece of Plasmodium, usually already attached to an oat flake, is 
placed in the capital city. The Petri dishes with Plasmodium were kept in darkness, at temperature 22-25 C°, 
except for observation and image recording (Fig. pi). Periodically (usually in 12 h or 24 h intervals) the dishes 



(a) Africa (b) USA 

Fig. 1. Experimental laboratory images of protoplasmic networks developed by slime mould P. polycephalum on maps. 



were scanned in Epson Perfection 4490. Examples of typical protoplasmic networks recorded in laboratory 
experiments are shown in Fig. [2] Detailed examples and scenarios of colonisation of various regions by the 
slime mould are provided in original papers [2]-[T4],[5T|. 

In our experiments we always inoculated Plasmodium of P. polycephalum in a capital of any particular 
country. This is because in the majority of cases — but not in all cases indeed — capitals are the most 
populated and industrially developed urban areas, and the road 'diffusion' in ancient times was usually 
originated from the capital city. For example, see a rough scheme of turnpikes in England (Fig. |3(a)| ), which 
demonstrates a classical growth pattern, typical for fungi, myxomecetes and bacterial colonies. Moreover, 
we can even see evidence of secondary growth from other cities, e.g. in Fig. |3(b)| we see the main turnpike 
network growing from London and several sub-networks growing from Bristol, Ross, Leominster, Worcester, 
and Manchester. 

Unlike bacterial colonies and fungal mycelia, the Physarum Plasmodium is able to shift the mass of 
its body plan by the adaptive assembly and disassembly of its protoplasmic tube network and utilise the 
transport of structural components within the network to the leading edge of growth. The morphology of 
the Plasmodium is thus less dependent on the initial inoculation site and the active zone of growth can move 
throughout the environment as the Plasmodium forages for food. 

Also, from our previous experimental studies we know that when Plasmodium is inoculated in every 
point of a given planar set, the protoplasmic network formed approximates the Delaunay triangulation of the 
set [1 48]. Neither of the motorway networks considered match the Delaunay triangulation of major urban 
areas, thus simultaneous inoculation in all urban areas would not bring any additional benefits. 

As with every living creature, the Plasmodium of P. polycephalum does not always repeat its foraging 
pattern. To generalise our experimental results we constructed a Physarum graph with weighted-edges. A 
Physarum graph is a tuple P = (U, E, w), where U is a set of urban areas, E is a set edges, and w : E — >• [0, 1] 
associates each edge of E with a frequency (or weight) of the edge occurrence in laboratory experiments. For 
every two regions a and b from U there is an edge connecting a and b if a Plasmodium's protoplasmic link is 
recorded at least in one of k experiments, and the edge (a, b) has a weight calculated as a ratio of experiments 
where protoplasmic link (a, b) occurred in the total number of experiments k. We do not take into account 
the exact configuration of the protoplasmic tubes but merely their existence. In original papers [2]-[T4],[5T| 




(d) Brazil 



(e) Canada 



(f) China 




(j) Mexico 



(k) The Netherlands 



(1) USA 



Fig. 2. Exemplar configurations of protoplasmic networks developed by slime mould P. polycephalum on major urban 
areas U obtained in experimental laboratory studies. 




Fig. 3. Turnpike road network, (a) A scheme of London turnpike network in England, c. 
(b) A scheme of the England turnpike network in 1750. Modified from [21] . 



1720. Modified from \T8\ 



we dealt with threshold Physarum graphs P(0) = (U, T(E), w,0). The threshold Physarum graph is obtained 
from the Physarum graph by the transformation: T(E) = {e G E : w(e) > 0}. That is, all edges with weights 
less than or equal to are removed. With the increase of in a family of threshold Physarum graphs {P(0), = 
0, 1, 2, • • • , k — 1} the graphs undergo the following transitions (see country-specific details in [2]-[T4].[5T]): 
non-planar connected — >> planar connected — >> disconnected — )> all nodes are isolated 

In the present paper we consider only 'stressed' Physarum graphs P(0'), which have maximum possible 
values 0' and yet remain connected: 0' = max{# : P(0) = connected}. These Physarum graphs are shown in 
Fig. pi Values of 0' for studied regions are illustrated in Fig. |6(a)| 

To compare the Physarum graphs with motorway graphs we construct a motorway graph H as follows. Let 
U be a set of urban regions, for any two regions a and b from U, the nodes a and b are connected by an edge 
(a, b) if there is a motorway starting in a and passing in the vicinity of b and not passing in the vicinity of any 
other urban area c G U. Motorway graphs extracted from maps of motorway /highway /expressway /autobahn 
networks are shown in Fig. [5] 

Motorway and Physarum graphs were compared directly and using integral measures and indices. Let m 
be a number of edges in motorway graph H and / be number of edges in Physarum graph P, and i and j 
be nodes, and M and F be adjacency matrices. Direct matching between motorway and Physarum graphs 
is calculated as \i = ^ J^-- £(M^ , Fij) , where £(Mij,Fij) = 1 if Mij = Fij, and 0, otherwise. An economy 
of matching is calculated as e = y. Also, we compared the graphs by their average shortest path measured 
in nodes, average shortest path measured in normalised edge lengths (for each edge e G E we normalised 
its Euclidean length /(e) as /(e) «— max su^\. e f e ^\ ^ average degrees (sum of degrees of nodes divided by a 
number of nodes), average edge length (of normalised edges), diameters (longest shortest path) in nodes and 
normalised edge lengths, and maximum number of vertex independent cycles (two cycles are independent of 
each other if they do not share nodes or edges) . 

To measure 'compactness' of graphs we calculated average cohesion: let d be an average degree of a graph 
G and v^ be a number of common neighbours of nodes i and j, and d{ is a degree of node z, then cohesion 
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Fig. 4. Physarum graphs P(0') for highest values of which do not make the graphs disconnected. 
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Fig. 5. Motorway graphs H. 
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Fig. 6. (a) 0' versus number of nodes, (d) Matching \i versus economy e. 



Kij between nodes i and j is calculated as K{j - 
index [42], 77-index [24], and Randic index [44] . 
The Harary index is calculated as follows: 77 



di~\-dn 



Three topological indices were calculated: Harary 



\ ^2ij x(Dij), where D is a graph distance matrix, where 
Dij is a length of a shortest path (in normalised edge lengths) between i and j, and x(^ij) — D^ 1 ifij^j and 
0, otherwise. The 77-index shows a relationship between the total length of the graph L(G) and the distance 



along its diameter D(d) [24]: 77 



D{d)' 



The Randic index [44] is calculated as R = ^- • Cij * ( 



where C^ is an adjacency matrix, C = M or C = F. 



■ yj(di*dj) J 



3 Results 

Matching and Economy Top three regions in the best matches \i between motorway and Physarum 
graphs are Malaysia, Italy and Canada and top three most economically e matched are Italy, Brazil and UK 

(Fig. log. 

Finding 1 Let C\ <C^ if \i(C\) < //(C^), then regions can be arranged in the following hierarchy of absolute 
Physarum matching: USA <\ Brazil <\ {Mexico, Iberia} <\ Australia <\ {Germany, UK } <\ {Africa, the 
Netherlands } <\ {China, Belgium } <\ Canada <\ Italy <\ Malaysia. 

We can consider a product of matching to economy uo = /i • e as a rough parameter for estimating 
'slime-optimality' of motorways approximation. By values of uo regions can be arranged in the descending 
slime-optimality as follow (exact values of uo are in brackets): 

1. Italy (0.85) 

2. Malaysia (0.76) 

3. Canada (0.67), China (0.65), Belgium (0.6), Africa (0.6), UK (0.59) 

4. Netherlands (0.49), Germany (0.48), Australia (0.47), Brazil (0.44), Mexico (0.42), USA (0.4) 

5. Iberia (0.34) 



Average degrees Averages degrees of motorway graphs are usually lower than degrees of the corresponding 



Physarum graphs (Fig. 7(a) ). This is particularly visible with USA and Brazil motorway graphs, which average 
degrees are nearly 4.5 while Physarum graphs approximating the motorways have twice less average degree. 
Belgium, Canada and Malaysia show almost perfect match between Physarum and motorway graphs in 
average degrees (Fig. |7(a~). 
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(f) Diameter in nodes 
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Fig. 7. (a) Average degrees, (b) Maximum number of vertex independent cycles, (c) Average edge link, (d) Average 
length of a shortest path in nodes, (e) Average length of a shortest path in real values, (f) Diameter in nodes. 
(g) Diameter in real values. 



Maximum number of independent cycles Motorway and Physarum graph of China, the Netherlands, 
Canada and Italy have the maximum number of independent cycles (Fig. |7(b)[ ) and show a good match. 
Other regions can be subdivided on two groups: 

— number of independent cycles is higher in motorway graphs: Africa, Australia, Brazil, Mexico, UK, USA 

— number of independent cycles is higher in Physarum graphs: Belgium, Berlin, Canada, Germany, Iberia 
and Malaysia. 

The maximum number of independent cycles may characterise two properties of transport networks: fault- 
tolerance (more cycles indicate more chances for a transported objects to avoid faulty impassable links, sites 
of accidents, and jams) and locality (distant links increase chances of two cycles of sharing a node). Thus, 
we can propose that China motorway network is most fault-tolerant and locally connected, while transport 
networks in Canada, Italy, Malaysia and UK could be sensitive to disasters and overloads. 

Average edge length The closest match between Physarum and motorway graphs in average edge length 
is observed for Italy, Iberia, Germany, Mexico, Canada and Malaysia (Fig. |7(c)| . Highest mismatch is shown 
by the Netherlands (edges of motorway graph are longer than edges of Physarum graph), and Brazil, Africa, 
and USA (edges of Physarum are longer, in average, than of motorway graphs). 

Average shortest paths The average length of a shortest path between nodes shows little match between 
Physarum and motorway graphs (Figs. [7(d)] and [7(e)] ) . Usually Physarum graphs exhibit 1.5-2 times longer 
average shortest paths, this varies however between regions. Malaysia and Africa are the regions with longest 



average shortest paths, measured in nodes, in motorway and Physarum graphs (Figs. 7(d)). When shortest 
paths are measured in normalised edge length, Africa and Brazil have longest average shortest paths in 
Physarum graphs, and Malaysia and Italy in motorways graphs. Countries which show closest match between 
Physarum and motorway graphs — in average shortest path, measured in nodes (Figs. |7(d)[ ) — are Canada, 



China, the Netherlands; and — in average shortest path measured in normalised lengths (Fig. 7(e) ) — are 
Belgium, Canada, China. 

Diameters Being the longest shortest path, the diameter shows even less matching between Physarum 



and motorway graphs (Figs. 7(f) and 7(g) ) than average shortest path matching. Physarum graphs match 
motorways in diameter measured in nodes for Belgium, and in diameter measured in normalised lengths for 
Belgium, China, and the Netherlands. 

Cohesion For most regions considered, average cohesion of Physarum graphs is typically higher than of 
motorway graphs (Fig. |8(a) ). The difference is particularly strong for Brazil and USA, e.g. average cohesion 



of Brazilian motorway network is 0.8 while of the corresponding Physarum graphs it is 0.4. There is a match 
between Physarum and motorway graphs of Canada, Italy, and the Netherlands. The top three entries with 
highest cohesion of motorway graph are Brazil, Canada, China, and the Netherlands, and the top three 
regions with highest cohesion of Physarum graphs are Brazil, China and USA. Cohes ion of each edge in the 
complete graph is 2 ?n-i) ~ 0-44, this limit is nearly approached by H (Brazil) (Fig. 



a)). Motorway and, 



especially, Physarum, graphs of Malaysia and Italy show minimal cohesion, because cohesion is zero in chains 
and graphs with cycles over three nodes. 

Finding 2 Physarum matches motorway network of Canada, Italy and the Netherlands in terms of com- 
pactness, local densities and fault-tolerance of transport networks. 

Average cohesion is an indicator of compactness [28 . A sub-set of the graph with high cohesion remains 
connected even when some edges are removed, thus cohesion may characterises stability [47 or even fault- 
tolerance of graphs. The cohesion of a node is the minimum number of edges whose deletion makes the node a 
cut node of the resulting graph [45 , thus the cohesion is used to characterise a local density of sub-graphs [55] , 
and it is related to centrality [20] , and statistical properties of connectivity of graphs [52] . 
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Fig. 8. (c) Average cohesion, (c) the Harary index, (d) the 77 index, (e) the Randic index. 



The Harary index The Harary index [42] is well know for its predictive properties in chemistry [29 30 , and 
the index is based on 'the chemists' intuitive expectation that distant sites in a structure should influence each 
other less than the near site' [37] • This is probably not the case with slime mould and man-made motorway 
networks. Only four out fourteen regions satisfy the relation |1 — ^f#y| < 0.1: China, Canada, Belgium and 
Malaysia, where rj is the Harary index. Physarum poorly approximates motorways, in terms of the Harary 
index, in Brazil and USA (Fig. |8(b)[ ). 



The 17-index Five regions show over 0.1 mismatch, |1 



n(P) 
n(H)\ 



> 0.1, between the 77- indices of Physarum 
)J. This result is quite 



and motorway graphs: Germany, Iberia, the Netherlands, Malaysia, and USA (Fig. 

interesting. Recall, that the 77-index is a ratio of a total length of all normalised edges of a graph to a 
distance along the graph's diameter. Physarum graphs neither match motorways in diameters (Figs. mfJ and 
7(g)| ), nor do we witness a good match in average edge length, or shortest paths (Figs. 7(c)[ |7(d)| and |7(e)| ). 
However, when these factors are considered in proportion the match between the graphs occurs. 



The Randic index 

Finding 3 Physarum perfectly approximates motorway networks in terms of the Randic index. 

The Randic index R shows impeccable match between Physarum and motorway graphs (Fig. |8(d)| and 
Tab. ll. The largest value 1 — 4pj = 0.07 is for Belgium motorway graph, and corresponding Physarum 

graph. Star-graph has a minimum the Randic index y/n — 1 [19] • For an arbitrary graph G the boundaries 
are \Jn — 1 < R(G) < |. As we can see in Tab. u\ indices for all motorway and Physarum graphs are very 
close to upper boundary. The highest values for motorways are in Italy and USA, and for Physarum graphs 
are in Belgium, Malaysia and Canada. 

The Randic index R (originally called by Milan Randic as molecular branching index) [44] characterises 
relationships between structure, property and activity of molecular components [31 . The index relates to 



Table 1. Values of the Randic index. 



Country 


R(H) R(P) 


fl(H) 

H(p) 


n 
2 


Africa 


16.98 16.88 


-0.006 


17.5 


Australia 


11.9 12.17 


0.022 


12.5 


Belgium 


9.53 10.27 


0.071 


10.5 


Brazil 


10.04 10.14 


0.01 


10.5 


Canada 


7.25 7.76 


0.066 


8 


China 


15.08 15.1 


0.002 


15.5 


Germany 


9.94 10.15 


0.02 


10.5 


Iberia 


10.62 11.11 


0.045 


11.5 


Italy 


5.22 5.22 





5.5 


Malaysia 


9.63 9.88 


0.025 


10 


Mexico 


8.87 8.95 


0.009 


9.5 


Netherlands 


10.09 10.17 


0.008 


10.5 


UK 


4.61 4.7 


0.021 


5 


USA 


9.71 9.41 


-0.031 


10 



diameter [25] and is actually the upper boundary of diameter [56]. It also relates to chromatic numbers of 
graphs and eigenvalues of adjacency matrices [36]. There are proven linear relations between the Randic index 
and molecular polarisability, cavity surface areas calculated for water solubility of alcohols and hydrocarbons, 
biological potencies of anaesthetics [34 , water solubility and boiling point [26] and even bio-concentration 
factor of hazardous chemicals [46]. Estrada [32] suggested the following structural interpretation: the Randic 
index is proportional to an area of molecular accessibility, i.e. area 'exposed' to outside environment. Or we 
can say that the index is inversely proportional to areas of overlapping between spheres of specified radius 
enclosing the nodes. The more overlapping the less is Randic index. In terms of transport networks, we can 
interpret external accessibility as transport inaccessibility, proportional to areas of country not served by 
existing motorway links. 

Finding 4 Physarum well approximates motorway graphs in terms of transport accessibility. 

Along the above discourse we can speculate that UK, Italy and Canada (first three regions with the 
smallest Randic indices) have better transport coverage of their territories than Africa, China and Australia 
(top three regions with the highest Randic indices). 



Table 2. Geographical regions with extremal values of measures over motorway graphs H and Physarum graphs P. 



Measure \i 


maxjLi(H) 


min/i(H) 


max/x(P] 


min/i(P) 


Average degree 


China 


Malaysia 


China 


UK 


Number of independent cycles 


China 


Malaysia 


China 


UK, Mexico 


Average cohesion 


Brazil 


Malaysia 


Brazil 


Malaysia 


Average edge length 


the Netherlands Malaysia 


Italy 


China 


Shortest path, nodes 


Malaysia 


UK 


Malaysia 


the Netherlands 


Shortest path 


Malaysia 


UK, USA 


Africa 


China 


Harary index 


China 


Italy 


China 


Italy 


77-index 


Africa 


UK 


Africa 


UK 


Randic index 


Africa 


UK 


Africa 


UK 


Diameter, nodes 


Malaysia 


Iberia, UK, 


USA Malaysia 


UK 


Diameter 


Malaysia 


USA 


Africa 


China 



Extremal regions We call a region extremal if it displays minimum or maximum values of at least one 
measure over its motorway or Physarum graphs. The extremal regions are listed in Tab. [2] 



Africa shows maximum 77-index and Randic index on both H an P, and maximum average shortest path 
and diameter on P. This might indicate on critical dependencies between geographically close urban 
areas, large territorial spread of transport networks and relatively higher density of urban area along 
coasts (Fig. |5(a)1 ). 

Brazil shows maximum average cohesion on H and P. This is because Brazil has highest (amongst regions 
studied) number of locally connected sub-graphs with larges number of dependent cycles (Fig. |5(d")| ). 
China shows maximum average degree, number of independent cycles and Harary index on H and P; 
minimum average edge length, average shortest path and diameter on P. These indicate on high accessi- 
bility of major urban area in China, and fault-tolerance of Chinese motorways at a large-scale (Fig. |5(f)| ). 
The expressway network in China has been developed much recently, known as the national trunk high- 
way system, and it is a high-standard transport system planned by the central government. The system 
is designed to be optimal and many factors were properly taken into account including terrains and 
landscapes. 

Iberia has minimum diameter in nodes on H. The man-made transport network structure resembles a 
wheel with an 'axle' at Madrid, most major urban areas around coast forming a 'rim' linked to Madrid 
by 'spokes' (Fig. |5(h)[ ). 

Italy shows minimum Harary index on H and P, and maximum average edge length on P. These are 
due to tree-like structure of the transport networks and constraining of the urban areas in the prolonged 
shape of the country (Fig. |5(i)| ). 

Malaysia shows maximum shortest path in nodes and diameter in nodes on H and P; minimum average 
cohesion on H and P; maximum average shortest path and diameter on H; minimum average degree, 
number of independent cycles, average edge length on H. This is because Malaysian transport network 
does not have cycles and consists of a chain connecting urban areas along western coast (Fig. |5(j)| . 
Mexico shows minimum number of independent cycles on P because the slime mould approximation is a 
tree, almost a chain with few branches (Fig. |4(k)| ). 

The Netherlands shows maximum average edge length on H and minimum average shortest path in 
nodes on P, due to relatively compact location of urban areas with high density of local transport links 
(Figs. [5(1)1 and [4(1)1). 

UK shows minimum 77-index and Randic index and diameter in nodes on H and P; minimum average 
shortest path in nodes, shortest path on H, and minimum average degree and number of independent 
cycles on P (Figs. 5(m) and 4(m)|); 



— USA shows minimum average shortest path, diameter in nodes and diameter on H. These are because 
the motorway system in USA was built with optimality yet efficiency in mind. 

Bio-rationality of measures In Tab. [3] we provide binary evaluation M(C, /i) of matching between 
Physarum P and motorway H graphs calculated for each country C and measure \i as follows: M(C, \i) = 1 
if 1 1 — f-p( C ll | < 0.1 and 0, otherwise. A bio-rationality /3 of a measure fi is a number of regions C for which 
/i(H(C)) = /i(P(C)) (Tab. [3} bottom row). 

Finding 5 Hierarchy of bio -rationality of measures: Randic index >p H >p Average edge length, Harary 
index >p Average degree, Number of independent cycles >p shortest path, diameter >p average cohesion >p 
diameter in nodes >p shortest path in nodes. 

Matching between motorway graphs and Physarum graphs is most strongly expressed in the Randic index; 
further measures amongst the top ones are 77-index, Harary index, and edge length. Thus we can enhance 
the Finding [3] as follows. 

Finding 6 The Randic index is the most bio- compatible measure of transport networks. 

The Randic index is used to characterise relations between structure, property and activity of chemical 
molecules [31] , thus we can speculate that in terms of structure-property-activity Physarum almost perfectly 
approximates motorway networks in all regions! 



Table 3. Matching M between Physarum and motorway graphs for each country and each measure. 





Av 

degree 


Cycles 


Co- 
hesion 


Edge 
length 


SP in 
nodes 


SP 


Harary- 
index 


Pi- 
index 


the 

Randic 

index 


Diame- 
ter, 
nodes 


Diame- 
ter 


Bio- 
ratio- 
nality 
of 

motor- 
ways 


Africa 

Australia 

Belgium 

Brazil 

Canada 

China 

Germany 

Iberia 

Italy 

Malaysia 

Mexico 

The Netherlands 

UK 

USA 


1 
1 

1 
1 


1 
1 

1 

1 


1 
1 


1 

1 
1 
1 
1 

1 




1 

1 
1 


1 

1 
1 

1 
1 

1 


1 
1 
1 
1 
1 
1 

1 

1 
1 


1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 


1 


1 
1 

1 


1 
2 
7 
3 
7 
6 
3 
2 
5 
5 
3 
4 
3 
1 


Bio-rationality of 
measures 


4 


4 


2 


6 





3 


6 


9 


14 


1 


3 





Bio-rationality of motorways We calculate a bio-rationality of a motorway graph as follows p(C) = 
^n l£(|l — (P(C)) I — e )' €(P) = 1 ^ predicate p is true, and otherwise. We chose e = 0.1. Values of 
bio-rationality of motorways are shown in last column in Tab. [3j 

Finding 7 Hierarchy of bio -rationality of regions is as follows: 
{Belgium, Canada} > p China > p { Italy, Malaysia } > p the Netherlands 
> p { Brazil, Germany, Mexico, UK } > p { Africa, USA } 

Comparing hierarchies of absolute matching, see Finding [I] with the above hierarchy of bio-rationality 
we find Belgium, Canada, China and Italy are at the intersection of first three levels of the hierarchies. We 
omit Italy because its shape is intrinsically restrictive and invokes rather trivial architectures of protoplasmic 
networks. 

Finding 8 Motorway networks in Belgium, Canada and China are most affine to protoplasmic networks of 
slime mould P. polycephalum. 



4 Discussion 



Based on the results of our previous laboratory experiments with slime mould imitating the development 
of transport networks in fourteen regions [SJ-IHLEI] we undertook a comparative analysis of the motorway 
and protoplasmic networks. We found that in terms of absolute matching between slime mould networks 
and motorway networks the regions studied can be arranged in the following order of decreasing matching: 
Malaysia, Italy, Canada, Belgium, China, Africa, the Netherlands, Germany, UK, Australia, Iberia, Mexico, 
Brazil, USA. We compared the Physarum and the motorway graphs using such measures as average and 
longest shortest paths, average degrees, number of independent cycles, the Harary index, the il-index and 
the Randic index. We found that in terms of these measures motorway networks in Belgium, Canada and 
China are most affine to protoplasmic networks of slime mould P. polycephalum. With regards to measures 
and topological indices we demonstrated that the Randic index could be considered as most bio-compatible 
measure of transport networks, because it matches incredibly well the slime mould and man-made transport 
networks, yet efficiently discriminates between transport networks of different regions. 
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